A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

06 (08.02.4, 08.05.1, 08.06.3, 10.11.1, 10.15.1, 10.15.2 (Bochum 2)) 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



OO 
ON 

o 

Q 



> ■ 

m : 
m . 

(N ■ 
(N ■ 

oo : 

ON ■ 
Or 

i ' 
o ■ 

a: 



Kinematics and binaries in young stellar aggregates. 
I. The trapezium system BD+00°1617 in Bochum 2* 

Ulisse Munari 1 and Lina Tomasella 2 

1 Osservatorio Astronomico di Padova, Sede di Asiago, 1-36012 Asiago (VI), Italy (munari@pd. astro . it) 

Osservatorio Astrofisico del Dip. di Astronomia, Universita di Padova, 1-36012 Asiago (VI), Italy (tomasella@pd. astro . it) 

Received date ; accepted date 



Abstract. Internal kinematics, spectroscopic binaries 
and galactic motion are investigated for the trapezium 
system BD+00°1617 (which lies at the heart of the young 
open cluster Bochum 2) by means on 73 high resolution 
Echelle+CCD spectra secured over the period 1994-98. 
Two of the three O-type member stars are found to be bi- 
naries on close and highly eccentric orbits of 6.8 and 11.0 
day period. The spectra of the two binaries show large 
variations in the half-intensity and equivalent widths of 
the Hel absorption lines, which are intrinsic to the pri- 
maries and are not correlated to the orbital phase. Astro- 
metric and radial velocities exclude that one of the com- 
ponent is leaving as a runaway star, but upper limits are 
still compatible with the trapezium evaporating at very 
low relative velocities. The projected rotational velocities 
of the three constituent O-type stars are low. This con- 
forms to expectations from the high frequency of binaries. 
The observed radial velocity of Bochum 2 agrees with the 
Hron (1987) expression for the Galaxy rotation inside the 
latter's quoted errors. 

Key words: Binaries: spectroscopic - Stars: early type 
- Stars: fundamental parameters - Galaxy: kinematics 
and dynamics - Open clusters and associations: individual 
(Bochum 2) 



1. Introduction 

With this paper we begin a series devoted to the results 
of a prolonged high resolution spectroscopic monitoring 
of early type members of young open clusters, trapezium 
systems and OB associations. The emphasis is on accu- 
rate and extensive radial velocities (secured over several 

Send offprint requests to: U. Munari 

* Tables 2 and 3 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) 
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years), suitable to trace (a) the internal kinematics, (b) 
the galactic motion, (c) the binary star content and the 
spectroscopic orbits. The high dispersion spectra are also 
useful to derive the projected rotational velocities, the 
chemical composition and the evolutionary status of the 
member stars as well as to study the interstellar medium 
inside and towards such young stellar aggregates. These 
latter aspects are however broadly outside the scope of the 
present series and will be considered elsewhere. 

Our program encompasses a selection of young stel- 
lar aggregates spanning a range in total mass, number 
of member stars, age, total energy, etc. The targets have 
been chosen, whenever possible, among those with already 
published high quality astrometry and photometry. We 
believe that much can be learned by extensive and homo- 
geneous high-res spectral campaign data on such targets 
and to this aim we begun in 1990 a long term program 
with the telescopes and spectrographs of the Asiago As- 
trophysical Observatory. The ambitious goal is to search 
for observationally proved grand-relations governing the 
dynamical evolution of the young stellar aggregates which 
could then be compared with the predictions of N-body 
simulations. Such grand-relations are expected to become 
apparent when a significant number of young stellar ag- 
gregates will have been studied in the most homogeneous 
way and for such a reason we will postpone later in the 
present series a throughout discussion of the results accu- 
mulated on the way. 

We start with a so far poorly studied trapezium system 
(BD+00°1617) at the heart of a dim and very young open 
cluster (Bochum 2), which lies at a great distance in the 
general anti-center galactic direction (at ~ 14 kpc galacto- 
centric distance). 

Trapezium systems are groups of three or more stars 
whose apparent separations onto the sky are similar, as 
for the prototype of the class, the four O stars at the core 
of the Orion open cluster and nebula. OB associations 
contain typically three to five trapezia and these con- 
figurations have a high frequency among T-associations 
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too. Within stellar associations, trapezia appear to clus- 
ter close to the nuclei of open clusters and emission neb- 
ulae. A trapezium-like configuration is in general dynam- 
ically unstable and therefore it must be young. It has to 
evolve either into (a) hierarchical systems (i.e. ones in 
which the separation of successive binaries increases by 
large factors; they are dynamically stable and are present 
among systems of all ages), or into (b) hard binaries (with 
possible ejection from the system of one or more of the 
trapezium components). Ambartsumian (1954) computed 
in lv2x 10 4 AU the largest stable dimension of a trapez- 
ium in the solar vicinity against perturbation by passing 
stars: the one order of magnitude larger extension (lxlO 5 
AU) for cataloged trapezia reinforces the notion that they 
must be dynamically young systems. The first catalog of 
trapezium systems has been compiled by Ambartsumian 
(1954), listing 108 entries and a more recent one (442 sys- 
tems) has been presented by Salukvadze (1978). Trapezia 
have been traditionally found by filtering catalogs of vi- 
sual binaries, like the Index Catalog of Visual Binary Stars 
(IDS) by Jeffers et al. (1963). 

Bochum 2 (a 20 oo = 06 h 49 m , S 2 q q = +00°23', I = 
212.3, b = -0.4) has been studied by Moffat & Vogt 
(1975), Turbide & Moffat (1993) and Munari & Carraro 
(1995, hereafter MC95). At its center lies BD+00°1617A- 
B-C, a system of three O-type stars which luminosity 
largely dominates the apparently dim surrounding clus- 
ter. The three O stars are plotted in the finding chart of 
Figure 1. They appear equally spaced and on a straight 
line in the projection onto the sky. 

Bochum 2 is not included in catalogs of trapezium sys- 
tems most probably because BD+00°1617A-B-C is oddly 
missing from catalogs of visual binaries. However, it fully 
satisfies all classification criteria and in the following it 
will be discussed as a validated trapezium system. The 
three components will be identified in the following ac- 
cording to the numbering by MC95 (#20, #21 and #23, 
cf. Figure 1). 

2. Observations 

The program stars are listed in Table 1. Photometric 
and spectroscopic data come from MC95. Positions are 
from the Astrographic Catalog for the epoch 1901.0 and 
equinox J2000.0 (see section 2.2), placed on the reference 
system of the Hipparcos Catalogue. The journal of the 
observations is given in Table 2. 
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Fig. 1. Finding chart for the program trapezium stars at the 
heart of the open cluster Bochum 2. Orientation and scale are 
given. Capital letters relates to BD+00°1617 sub-classification, 
the numbers to the identification in Munari & Carraro (1995). 
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2.1. Spectroscopy 

Seventy-three high resolution spectra of the program stars 
have been secured from Nov. 1994 to Feb. 1998 with the 
1.82 m telescope and Echelle+CCD spectrograph of the 
Astronomical Observatory of Padova at Asiago (Cima 
Ekar). The wavelength range covered on a single frame 
was generally from RS to Hcl 6678 A, with some spectra 



Fig. 2. Change with time of the distance between stars A and 
B compared to that between stars B and C [R = A(AB / BC)]. 
The dashed line is a last square fit and the error bars are 3<7 
values. Solid circles are from the Astrographic Catalogue, the 
USNO Twin Astrographic Catalogue and some recent CCD 
observations, the open circles from Palomar Schmidt POSS-I 
(USNO A1.0 Catalogue) and GSC-I. 
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Table 1. Program stars. Identification numbers, photometry and spectroscopy from Munari & Carraro (1995; MC95), positions 
from the Astrographic Catalog 2000 reduced to the Hipparcos Catalogue reference system. 



MC95 


BD 


LS 


CSC 


V 


U-B 


B-V 


spectrum 




QJ2000 <5.72()()() 


# 






# 












(epoch 1901.0) 


20 


+00°1617 C 


VI +00 26 


1481994 


11.21 


-0.52 


0.49 


09 V 


06 h 


48 m 51. s 300 +00° 22' 21. "49 


21 


+00°1617 B 


VI +00 25 


1482125 


11.01 


-0.48 


0.53 


07 V 


06 h 


48 m 50. s 484 +00° 22' 37."62 


23 


+00° 1617 A 


VI +00 24 


1481965 


11.38 


-0.52 


0.51 


09 V 


06 h 


48 m 49. s 541 +00° 22' 52. "71 




Fig. 3. Orbital solution for star #20. The solid line refers to 
the orbital solution given in Table 4. The dashed line is the 
baricentric velocity (69.8 km sec -1 ). 



Fig. 4. Orbital solution for star #21. The solid line refers to 
the orbital solution given in Table 4. The dashed line is the 
baricentric velocity (63.0 km sec -1 ). 



extending up to 9900 A and a few others down to 3865 
A (cf. Table 2 for details). The instrumental PSF (from 
the FWHM of night-sky lines and thorium lines in the 
comparison spectra) is generally uniform over the various 
order and pretty constant during the 5 years of observa- 
tions, corresponding to an average of 14 km sec -1 for full 
resolution frames and to 20 km sec -1 for those with a 
2x2 binning (corresponding respectively to 0.22 and 0.32 
A resolution at H/3). Details on the instrument, its perfor- 
mances and the use are given in Munari & Zwitter (1994). 

The spectra have been extracted and calibrated in a 
standard fashion using the IRAF reduction package. Great 
care has been taken to ensure the highest quality in the 
wavelength calibration and therefore its suitability for ac- 
curate radial velocity measurements. The flexure pattern 
of the Asiago spectrographs (and their negligible influ- 
ence on derivation of accurate RVs) has been studied in 
detail by Munari & Lattanzi (1992). The small (however 
always smooth and predictable) flexures induced by tele- 



scope tracking on the targets has been removed accord- 
ingly. 

Some checks have been performed to control the accu- 
racy of the wavelength scale of the calibrated 1-D spectra. 
A brief description of the results from these checks follows. 

(i) At the S/N ratios given in Table 2 only a few of 
the strongest telluric absorption lines are detectable on 
a single spectrum. Therefore, for each program star we 
have summed up in the velocity space all wavelength cali- 
brated 1-D spectra. On the summed-up spectra many tel- 
luric lines (as weak as to produce an equivalent width of 
only 0.004 A) are measurable. Their mean radial velocity 
has turned out to be 0.0 km sec -1 with a dispersion of 0.52 
km sec -1 , largely dominated by the difficulty to measure 
so weak features. This supports the conclusion by Munari 
& Lattanzi (1992) that the Asiago Echelle+CCD spectro- 
graph has a completely negligible (if any) spectrograph 
velocity or systematic offset. 
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Table 2. Journal of observations. N = spectrum number from the telescope logbook; # = program star identification number; 
HJD — heliocentric JD - 2400000; Res: f = full-resolution spectrum, b — 2 x 2 binned spectrum; S/N — signal-to-noise ratio 
for the stellar continuum close to the Nal D doublet (center of Echelle order 38), or center of order 26 for last entries. Column 
Orders lists the Echelle spectral orders covered on the given spectrum, which corresponds to the following wavelength ranges: 
23-29 = 9900-7600 A, 28-38 = 8130-5800 A, 29-41 = 7850-5375 A, 33-51 = 6900-4320 A, 33-52 = 6900-4240 A, 34-51 = 
6695-4320 A, 34-50 = 6695-4410 A, 34-53 = 6695-4160 A, 37-50 = 6150-4410 A, 38-57 = 5990-3865 A, 39-58 = 5835-3800 
A. The HJD of the last three spectra is the same because their spectra have been recorded simultaneously on the same CCD 
frame. 



N 


# 


Date 


HJD 


Res. 


Orders 


S/N 


17522 


21 


13.11.94 


49669.554 


f 


38- 


-57 


29 


17524 


20 


13.11.94 


49669.584 


f 


38- 


-01 


26 


17526 


23 


13.11.94 


49669.622 


f 


38- 


-57 


27 


17869 


21 


16.12.94 


49703.411 


f 


on 

IM- 


A 1 

41 


33 


17871 


20 


16.12.94 


49703.438 


f 


29- 


A 1 

41 


35 


17873 


23 


16.12.94 


49703.472 


f 


29- 


-41 


32 


18336 


21 


15.01.95 


49733.379 


f 


O A 

34- 


-53 


39 


18337 


20 


15.01.95 


49733.410 


f 


O A 

34- 


53 


48 


18339 


23 


15.01.95 


49733.457 


f 


34- 


-53 


47 


21191 


23 


10.12.95 


50061.650 


b 


O A 

34- 


-53 


33 


21237 


23 


11.12.95 


50062.589 


b 


O A 

34- 


-53 


38 


21239 


21 


11.12.95 


50062.602 


b 


34- 


-53 


39 


21241 


20 


11.12.95 


50062.616 


b 


O A 

34 


-53 


29 


21276 


23 


12.12.95 


50063.558 


b 


34- 


-53 


32 


21595 


23 


31.03.96 


50174.372 


b 


O A 

34- 


-53 


18 


21597 


21 


31.03.96 


50174.389 


b 


Q A 

34- 


-53 


15 


21696 


20 


06.04.96 


50180.314 


b 


34- 


-53 


29 


21698 


20 


06.04.96 


50180.334 


b 


34- 


-53 


23 


21700 


21 


06.04.96 


50180.368 


b 


34- 


-53 


26 


22588 


23 


30.12.96 


50447.578 


b 


33 


-52 


41 


22590 


21 


30.12.96 


50447.591 


b 


33- 


-52 


57 


22592 


20 


30.12.96 


50447.607 


b 


33- 


-52 


35 


23061 


20 


27.01.97 


50475.501 


b 


34- 


-53 


30 


23063 


21 


27.01.97 


50475.514 


b 


34- 


-53 


34 


23065 


23 


27.01.97 


50475.527 


b 


34- 


-53 


26 


23168 


20 


28.01.97 


50476.505 


b 


34- 


-53 


23 


23170 


21 


28.01.97 


50476.521 


b 


34- 


-51 


24 


23172 


23 


28.01.97 


50476.538 


b 


34- 


51 


19 


23644 


20 


17.02.97 


50497.358 


b 


33- 


51 


32 


23645 


21 


17.02.97 


50497.381 


b 


33- 


51 


36 


23647 


23 


17.02.97 


50497.404 


b 


33- 


-51 


36 


23723 


23 


18.02.97 


50498.353 


f 


33- 


-50 


23 


23725 


21 


18.02.97 


50498.401 


f 


33- 


-50 


16 


23823 


21 


21.02.97 


50501.342 


f 


34- 


-53 


38 


23825 


20 


21.02.97 


50501.390 


f 


34- 


-53 


36 


23827 


23 


21.02.97 


50501.437 


f 


34- 


-53 


29 


23899 


23 


22.02.97 


50502.330 


f 


34- 


-53 


40 



N 


# 


Date 


HJD 


Res. 


Orders 


S/N 


23901 


21 


22.02.97 


50502.365 


b 


34- 


-53 


49 


23903 


20 


22.02.97 


50502.389 


b 


34- 


-53 


50 


24457 


21 


25.03.97 


50533.367 


b 


34- 


-53 


21 


24459 


20 


25.03.97 


50533.384 


b 


34- 


-53 


18 


24461 


21 


25.03.97 


50533.400 


b 


34- 


-51 


17 


24477 


21 


27.03.97 


50535.266 


b 


34- 


-53 


32 


24479 


20 


27.03.97 


50535.278 


b 


34- 


-53 


35 


24483 


21 


27.03.97 


50535.306 


b 


34- 


-53 


21 


24485 


20 


27.03.97 


50535.322 


b 


34- 


-53 


32 


24487 


23 


27.03.97 


50535.346 


b 


34- 


-53 


36 


24489 


21 


27.03.97 


50535.369 


b 


34- 


-53 


28 


24491 


20 


27.03.97 


50535.386 


b 


34-53 


24 


24651 


20 


16.04.97 


50555.282 


b 


34- 


-50 


16 


24653 


21 


16.04.97 


50555.295 


b 


34- 


-53 


34 


24655 


20 


16.04.97 


50555.312 


b 


34- 


-53 


31 


24699 


20 


17.04.97 


50556.297 


b 


34- 


-53 


45 


24701 


21 


17.04.97 


50556.315 


b 


34- 


-53 


49 


24703 


23 


17.04.97 


50556.332 


b 


34 


-53 


48 


24770 


20 


18.04.97 


50557.290 


b 


34- 


-53 


39 


24772 


21 


18.04.97 


50557.313 


b 


34- 


-53 


42 


24774 


20 


18.04.97 


50557.332 


b 


34- 


-53 


26 


26468 


20 


14.11.97 


50766.603 


b 


34- 


-53 


38 


26472 


21 


14.11.97 


50766.650 


b 


34- 


-53 


45 


26474 


20 


14.11.97 


50766.673 


b 


34- 


-53 


42 


26506 


20 


15.11.97 


50767.595 


b 


34- 


-53 


30 


26510 


21 


15.11.97 


50767.637 


b 


34- 


-53 


32 


26512 


20 


15.11.97 


50767.661 


b 


34- 


-53 


25 


26518 


20 


15.11.97 


50767.700 


b 


34- 


-53 


25 


27068 


20 


13.12.97 


50795.522 


b 


28- 


-38 


24 


27070 


21 


13.12.97 


50795.546 


b 


28- 


38 


31 


27471 


20 


11.01.98 


50825.461 


b 


39- 


-58 


54 


27473 


21 


11.01.98 


50825.485 


b 


39- 


-58 


83 


27474 


20 


12.01.98 


50825.515 


b 


39- 


-58 


72 


27653 


20 


05.02.98 


50850.450 


b 


23- 


-29 


45 


27653 


21 


05.02.98 


50850.450 


b 


23- 


-29 


60 


27653 


23 


05.02.98 


50850.450 


b 


23- 


-29 


38 



(ii) Our spectra have been exposed for an average of 
30 minutes and at an average height over the horizon of 
38 degrees. In these condition strong night-sky lines are 
expected and have been indeed recorded on our spectra. 
On all spectra we have measured the three strongest night- 
sky lines, which lay on separate Echelle orders. The results 
confirm the high control over the wavelength calibration, 
with a mean radial velocity of 0.0 km sec -1 and a typical 



dispersion of 0.19 km sec -1 (the lower dispersion com- 
pared with the weak telluric lines has to be ascribed to 
the higher S/N ratio of the night-sky lines). 

(Hi) Finally, nearly all spectra cover the region of the 
interstellar Nal D doublet. We have measured its radial 
velocity, with the following results: 

J?V^ oI (#20) = 40.94 ±0.18 (a = 0.92) km sec _1 (l) 
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Table 3. Heliocentric radial velocities and associated errors for the three program stars. N°— spectrum number as in Table 2. 



#20 #21 #23 



N 


HJD 


RV^ 

XXV 




N 


HJD 


RV^ 

XXV 




N 


HJD 


RV^ 

XXV 




17524 


49669 584 


54.7 


1.3 


17522 


49669 554 

t: t7 \7 \7 17 . < J < J t: 


83.8 


1.2 


17526 


49669 622 


68.1 


2.9 


17871 


49703 438 

4:t7 1 U0.4:00 


56.0 


0.4 


17869 


49703 41 1 


54.7 


2.7 


17873 


49703 479 


68.3 


0.3 


18337 


49733 41 


44.0 


1.1 


18336 


49733 379 

Tit7 1 ' )• } • ' } 1 U 


133.7 


1.5 


18339 


49733 457 


69.2 


2.2 


21241 


50062 61 6 


48.9 


1.4 


21239 


50062 602 


71.4 


1.3 


21191 


50061 650 

< J V7 V7 V7 i . VJ O V J 


65.5 


3.6 


21696 


50180 314 


101.5 


2.9 


21597 


501 74 389 


135.2 


1.5 


21237 


50062 589 

'J V7VJ V7^ .UU<7 


67.0 


1.0 


21698 


50180 334 


113.4 


0.4 


21700 


50180 368 


7.3 


2.0 


21276 


50063 558 


68.5 


1.1 


22592 


50447 607 

• 7 V / 1 i J - V 7 ' I 


103.9 


1.7 


22590 


50447 591 


34.2 


2.5 


21595 


50174 372 


65.8 


3.6 


23061 


50475 501 


107.3 


0.8 


23063 


50475 514 


46.9 


0.7 


22588 


50447 578 


67.6 


1.9 


23168 


50476 505 


89.3 


6.0 


23170 


50476 521 


47.5 


5.4 


23065 


50475 527 


68.0 


3.0 


23644 


50497 358 


78.5 


0.9 


23645 


50497 381 


69.9 


0.7 


23172 


50476 538 

' 7 v 7 t I \l . 7 " 7 '. 7 


70.3 


5.5 


23825 


50501 390 


46.7 


1.4 


23725 


50498 401 


51.3 


1.1 


23647 


50497 404 


69.2 


1.3 


23903 


50509 38Q 


85.1 


1.8 


93893 


50501 349 


15.8 


1.9 


23723 


50498 353 


67.7 


2.5 


24459 


50533 384 


47.5 


0.9 


23901 


50502 365 

• 7 V ' • 7 V ' — . ' 7 V' ' 7 


22.1 


1.1 


23827 


50501 437 

' 7 V7 ' 7 V7 1 . i "7 1 


67.8 


2.2 


24479 


50535 278 


44.8 


0.8 


24457 


50533 367 


3.4 


3.6 


23899 


50502 330 


71.5 


2.7 


24485 


50535 322 


44.3 


3.2 


24461 


50533 400 


5.2 


3.8 


24487 


50535 346 


69.2 


1.0 


24491 


50535 386 


50.2 


1.2 


24477 


50535 266 


13.5 


3.5 


24703 


50556 332 


68.8 


2.1 


24651 


50555 282 


50.7 


0.4 


24483 


50535 306 


20.7 


2.5 


27553 


50850 450 


69.0 


7.0 


24655 


50555 31 2 


48.0 


0.8 


24489 


50535 369 


6.0 


3.6 










24699 


50556 297 


51.7 


4.9 


24653 


50555 295 


13.0 


2.2 










24770 


50557.290 


106.4 


0.8 


24701 


50556.315 


8.9 


5.0 










24774 


50557.332 


97.3 


0.4 


24772 


50557.313 


17.9 


1.1 










26468 


50766.603 


50.3 


6.0 


26472 


50766.650 


12.5 


3.9 










26474 


50766.673 


53.7 


2.9 


26510 


50767.637 


39.0 


4.6 










26506 


50767.595 


52.8 


1.1 


27070 


50795.546 


51.0 


2.4 










26512 


50767.661 


50.9 


4.7 


27473 


50825.485 


138.9 


2.0 










26518 


50767.700 


47.7 


3.1 


27553 


50850.450 


51.9 


5.1 










27068 


50795.522 


48.9 


7.4 


















27471 


50825.461 


115.5 


15.0 


















27474 


50825.515 


110.8 


1.5 


















27553 


50850.450 


49.9 


8.0 



















Table 4. Orbital elements for stars #20 and #21. The last row gives the weighted deviation of the observed radial velocities 
from the computed orbital solution. The quoted errors are the formal errors of the orbital solution. HJD = heliocentric JD - 
2400000. 



#20 #21 



period 


(days) 


6.858 


±0.004 


11.030 


±0.006 


baricentric velocity 


(km sec -1 ) 


69.8 


± 1.5 


63.0 


± 1.6 


semi-amplitude 


(km sec" 1 ) 


33.3 


± 2.3 


64.9 


± 2.6 


eccentricity 




0.43 


± 0.17 


0.21 


± 0.05 


asini 


(AU) 


0.0019 


± 0.0001 


0.065 


± 0.003 


To 


(HJD) 


50831.786 


± 0.278 


50835.749 


± 0.434 


Q 


(deg) 


316 


± 19 


322 


± 14 


mass function 




0.019 




0.29 




deviation 


(km sec" 1 ) 


4.2 




4.8 
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Table 5. Dispersion (a), mean and error of the mean for equivalent widths (EW) and half-intensity widths (HIW) of A 5876 A 
Hel line. The projected rotational velocities are from Equation (8). 



# 




EW (A) 


HIW (A) 


Vrot sini 




a 


mean 


a 


mean 


(km sec -1 ) 


20 


0.20 


0.73±0.04 


0.53 


3.01±0.13 


92±8 


21 


0.11 


0.60±0.03 


0.62 


3.14±0.14 


98±9 


23 


0.010 


0.706±0.003 


0.07 


1.54±0.02 


30±5 
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Fig. 5. Half-intensity width (A) of A 5876 A Hel on spectra Fig. 6. Equivalent width (A) of A 5876 A Hel on spectra of 
of programme stars ##20 (upper panel), ##21 (center) and ##23 programme stars ##20 (upper panel), ##21 (center) and ##23 



(lower panel). The dashed lines represent the mean values. 



(lower panel). The dashed lines represent the mean values. 



RV£ aI {#21) = 41.06 ±0.23 {a = 1.09) km sec -1 (2) 
RVg aI (#23) = 40.33 ±0.24 (a = 0.95) km sec -1 (3) 

where the standard error of the mean and the dispersion 
of the measurements are given. Again, the low values for a 
confirm the accuracy and repeatability of the wavelength 
scale during the night and over several years. 

2.2. Astrometry 

Some constraints on the internal dynamics of the trapez- 
ium system BD+00°1617 can be derived from data 
available in published astrometric catalogues, supple- 
mented by some recent ground-based CCD observations. 
BD+00°1617 was included in the first version of the Ty- 
cho Input Catalogue (TIC), but unfortunately it did not 
pass the Recognition Process after one year of satellite 
operation which led to the final version of the TIC (Halb- 
wachs et al. 1994). Plates from the POSS-II covering the 
Bochum 2 area have not yet been digitized by the time 



this paper has been written, and so preliminary GSC-II 
data are not available. 

The oldest epoch data for BD+00°1617 come from the 
Astrographic Catalogue as reduced to the ACRS released 
by the United States Naval Observatory (USNO, cf. Ur- 
ban & Corbin 1994, 1996). Bochum 2 is in the zone pho- 
tographed in 1901 by the Observatory of Algiers. Interme- 
diate epoch data are available from (a) the 1984 Palomar 
Schmidt plates used for the Guide Star Catalogue (GSC- 
I; Lasker et al. 1990), (b) the astrometric reduction of the 
plates of the original Palomar Sky Survey (POSS-I) per- 
formed at USNO and known as the A1.0 Catalogue and, 
finally, (c) positions measured along the USNO Twin As- 
trographic Catalogue project (Zacharias et al. 1996). 

Recent epoch astrometric data come from MC95 (the 
X,Y in their Table 3). Similar observations have been se- 
cured with the 1.22 m reflector of the Asiago Astrophysical 
Observatory during the commissioning in Feb. 1997 of a 
new CCD camera for the Newton focus (6 m focal length) . 
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For both MC95 and Asiago data, no field stars suitable to 
arrange a local astrometric reference system entered the 
CCD field around Bochum 2, and therefore no reduction 
to J2000.0 positions has been possible. Only relative po- 
sitions between components A, B and C of BD+00°1617 
have been measured. They are plotted in Figure 2 as the 
ratio of the distance from component A to B, with respect 
to distance from B to C. The same ratio has been com- 
puted for the other data sets above described and plotted 
in the figure. 

A 3<T error bar can be derived for the Astrographic Cat- 
alogue and the CCD data, but not unambiguously for the 
CSC and POSS ones. The latter two appear in Figure 2 
to be of poorer astrometric quality (possibly due to the 
blending with nearby faints stars of the associated clus- 
ter and the notoriously complicated focal field geometry 
of Schmidt telescopes) and will not be considered further 
on. 

Figure 2 does not support a detectable change of the 
relative distances of the three program stars (AB ~ BC ~ 
20. "5) over the last century. A least square solution (of 
quite low statistical significance) may indicate a change by 
~0.2% over the time span by the data in Figure 2 (about 
a century). At a distance (D) to Bochum 2 of <~ 6 kpc it 
corresponds to a projected relative velocity of 

V P rojected ~ 10 (jr^^j km sec" 1 (4) 

It does not support a runaway status for any of the pro- 
gram stars and it is largely compatible inside the er- 
rors with a negative total energy (a gravitationally bound 
trapezium). 

In view of the limited accuracy and overall paucity 
of available data, a confirmation and refinement of the 
astrometric results would be desiderable through a search 
in plates archives for suitable old and intermediate epoch 
astrometric plates covering the region of Bochum 2. 

3. Radial Velocities 

Radial velocities have been first measured on the cali- 
brated spectra by fitting individually each absorption line 
(with reference wavelengths taken from Moore 1959). This 
soon led to the recognition of star #23 as a constant radial 
velocity star and the other two program stars as binaries. 
The heliocentric radial velocity of star #23 (cf. Table 3) 
has been found to be: 

RV* 23 = 68.3 ± 0.3 (cr = 1.4) km sec" 1 (5) 

Three times the cr = 1.4 km sec -1 can be taken as the 
threshold for detection of binarity among the program 
stars. Using #23 as a RV standard star we have proceeded 
to re-evaluate radial velocities of #20 and #21 by cross- 
correlation (with the IRAF task fxcor). The results are 
listed in Table 3. 



Taking a mean between the velocity of #23 and the 
baricentric velocities of #20 and #21 (see next section and 
Table 4), the trapezium (and hence cluster) heliocentric 
radial velocity is (weighted mean): 

RVg = 68 ± 3 km sec" 1 (6) 

To the best of our knowledge, only Jackson et al. (1979) 
have so far published a radial velocity determination for 
Bochum 2. They reported a 49 km sec -1 from a few spec- 
tra of one or more of our same program stars. Unfortu- 
nately they gave no further detail. A ±7 km sec -1 er- 
ror is attached to this value by Turbide & Moffat (1993). 
Given the large amplitude of radial velocity variation by 
the the program stars and the absence of informations on 
the number and quality of the measurements by Jackson 
et al. (1979), their RV for Bochum 2 will not be considered 
in the following. 

4. Spectroscopic Orbits 

The large variability of the radial velocities of #20 and 
#21 compared to the constancy for #23 (cf. Table 3) ar- 
gues for #20 and #21 to be binaries. We have therefore 
proceeded with period search and determination of the 
spectroscopic orbits for the latter two stars. 

The orbital periods of #20 and #21 turned out to 
be 6 d 858 and ll d 030, respectively. The orbital solutions 
computed from these orbital periods are listed in Table 4 
and in graphical form are presented in Figures 3 and 4. 
The errors of the orbital solutions are of the order of 5% 
or less. The larger errors for the eccentricities may results 
from the uneven distribution of the observations along the 
orbital phase. 

5. Rotational Velocities 

To estimate the projected rotational velocities of the pro- 
gram stars (V rot sini) we have selected a sample of early 
type objects from the rotational velocity standards given 
by Sletteback et al. (1975). We have observed them with 
the same Echelle+CCD set-up used for BD+00°1617. 

Helium lines are notoriously a good choice for rota- 
tional velocity determination in early type stars (they are 
less affected by Stark broadening compared to hydrogen 
and are stronger than metallic lines). A 5876 A is the best 
recorded helium line in our spectra of BD+00°1617 and 
the half-intensity width (HIW) of the A 5876 A Hel line 
has been measured on all standard stars. In Figure 7 the 
HIW of the A 5876 A Hel line for the Sletteback et al. 
(1975) standard stars is plotted against the tabulated ro- 
tational velocity. The correction for the instrumental PSF 
has been applied to the data assuming that the PSF and 
the rotationally broaden profile (at low V ro t) combine as: 

a rotation = ° observed ~ °PSF CO 
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Fig. 7. Calibration for our Echelle+CCD spectrograph of the 
relation between half-intensity width (HIW) of the A 5876 A 
Hel line and rotational velocity (standard stars from Sletteback 
et al. 1975). The solid line is the least square fit given at the top 
of the figure (cf. Eq. 8). Spectral types and luminosity class of 
the observed stars are given next to their symbols. HIW errors 
do not exceed the symbol dimension. For errors in Vrot sinz see 
Sletteback et al. (1975). 



The instrumental PSF has been derived from the pro- 
file of night-sky emission lines (which follow the same 
optical path as the stellar light). The effect of the cor- 
rection is quite small and the induced displacement in 
Figure 7 is barely noticeable only for the two stars with 
V ro t sini < 50 km sec -1 . 

As expected, a quite straight linear relation exists in 
Figure 7 between the projected rotational velocity and the 
HIW, equally valid for main sequence and giant stars. A 
least square fit gives: 



V rot sim = 42.42 x HIW - 35 



km 



(8) 



The average HIW of the A 5876 A Hel line for the 
BD+00°1617 stars as measured on all available spectra 
is given in Table 5 together with the corresponding pro- 
jected rotational velocity inferred from relation (8). 

The low values of V rot sim for BD+00°1617 stars agree 
with the high fraction (2/3) of binaries, a general trend 
first noted by Abt & Hunter (1962 - and confirmed by 
later investigations, cf. Hack & Struve 1969) who pointed 
out how tidal energy dissipation in close binaries is ex- 
pected to act as a breaking mechanism on stellar rotation. 

6. Helium lines variability 

The equivalent width (EW) and half-intensity width 
(HIW) of the A 5876 A Hel line have been measured on all 



spectra of the program stars. The individual measures are 
plotted as function of time in Figures 5 and 6. In Table 5 
the dispersion of the measurements, their mean and the 
error of the mean are given. There is a striking evidence: 
Hel lines in the two binary stars vary by a large amount 
both in equivalent width and in broadness while they are 
remarkably constant in the non-binary #23 star. The am- 
plitude of variation is one order of magnitude the error 
of the measurements (fixed by the dispersion around the 
mean for the non-binary #23) and therefore the detection 
has to be considered as quite robust. 

The variation can in principle either be (a) intrinsic to 
the primary, or (b) ascribed to the visibility of the spec- 
trum of the secondary blended with the primary's one. 
The following discussion leads to the conclusion that the 
variations are most probably intrinsic to the primary. 



co 2 



CD 



x 



50 100 150 

radial velocity (km sec -1 ) 

Fig. 8. Half-intensity width (HIW) of Hel 5876 A plotted 
against radial velocity for #20 (filled circles) and #21 (crosses). 
The expected distribution for an unresolved double lined spec- 
troscopic binary is V-shaped, with maximum HIW reached at 
extrema of the velocity range. This is clearly not the case for 
both stars which distributions on the graph do not show any 
correlation with the radial velocity. 



First, we discuss the possibility that the variations are 
induced by the visibility of the spectrum of the secondary. 
Both the EW and HIW of the Hel line vary by a factor of 
two. Thus the intensity of the Hel line in the companion 
star should be similar to that of the primary and the two 
must have similar brightness. From Jasheck & Jaschck 
(1987) a general upper limit to the intrinsic rotational 
velocity of late O stars may be taken to be V rot ^300 km 
sec -1 . Comparison with the projected rotational velocity 
in Table 3 gives a lower limit to the inclination of the 
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rotational axis. The assumption that rotational axis and 
the orbital axis are roughly aligned puts the constraint 
on the mass of the companion from the mass function in 
Table 4. For star #20 this is 



,•#20 



#20 
orb 



> 18° 



Ml 



(Mr + M 2 ) 2 
M* 20 <8 M G 



< 0.64 



#20 

© 



(9) 



The masses are taken as M (07 V) = 28 M and M(09 V) 
= 19 M Q (cf. Conti 1975). This maximum mass for the sec- 
ondary in #20 correspond to a B5 main sequence, which 
is vastly too dim to rival in brightness and Hel line in- 
tensity with the 09 primary. The Hel line variability in 
#20 therefore cannot be ascribed to the visibility of the 
companion and it is therefore intrinsic to the primary. 

For #21 (an 07 V from Table 1) the situation remains 
unsettled because the same above line of reasoning on in- 
clination leads to 



.•#21 



i* 21 > 19° 

l orb =- ly 



Ml 



(Mr + M 2 )2 



M* 21 <31 Mi 



< 8.4 



#21 

© 



(10) 



which is compatible with a companion of the same mass 
and therefore the same brightness. 

The most stringent argument on the nature of the Hel 
line variability comes however from the search of a corre- 
lation between extrcma in the Hel variation and extrema 
in the radial velocities. The result is presented in Figure 8. 
Maximum split between the lines of the primary and sec- 
ondary (and therefore a maximum HIW value) has to be 
expected when the primary is passing through radial ve- 
locity extrema, resulting in a V-shaped distribution of the 
points in Figure 8. The latter clearly shows this not to be 
the case for both #21 and #20, which Hel line variability 
is intrinsic to the O primary star and appears triggered 
by the binarity. 

Absorption line variability is known to occur in early 
type stars, both single and binary systems, although no 
satisfactory explanation yet exists for the phenomenon. 
Jaschek, Jaschek and Kucewicz (1968) reported that in 
HD 125823 (B7 III) the Hel lines vary enormously from 
those of a B2 to those of a B9 star, while all other lines 
remain constant and photometric and RV variations are 
at most marginal (Norris 1971). We too (Lattanzi, Mas- 
sone and Munari 1992) identified two helium variable stars 
(their mean spectra suggesting a classification toward late 
main sequence B stars) in the young open cluster NGC 
225. Struve-Sahade effect is termed by Howarth et al. 
(1997) the systematic variation in relative line strength as 
a function of orbital phase in double-lined spectroscopic 
binaries of early types. In their investigation based on IUE 
spectra the frequency of incidence of such an effect on 
double-lined OB spectroscopic binaries is 1 in 4. When 
the effect is present, radial velocity curves appear more 



"noisy" than expected solely on the base of accuracy of 
the measurements (cf. Struve 1948 for the template case 
of the double-lined 08 V spectroscopic binary HD 47129, 
the Plaskett's star). 

7. Galactic rotation 

The component of the radial velocity due to the galactic 
rotation may be written as (cf. Hron 1987): 

RVq = w sin6 + u cosZ cos6 — v s'ml cos6 

- 2[A(R - R ) + a(R - R ) 2 } sinZ cos6 (11) 



where 



A : Oort' s constant = — 



R ( duj 



2 \dR 



t t Ro fd 2 uj\ 

a : curvature term = — — ■ 

4 \dR 2 ) Ro 



R = y/Rl + d 2 - 2R Q dcosl 

and R, R arc the cluster and Sun galacto-centric dis- 
tances, d the distance cluster-Sun, (l,b) the heliocentric 
galactic coordinates of the cluster and (u ,v ,w ) is the 
solar motion vector. Adopting from MC95 a cluster dis- 
tance of 6 kpc and (u ,v ,w ) = (-9.32,11.18,7.0) km 
sec -1 from Pont et al. (1994), Eq.(ll) rewrites as: 



RVq 



13.8 + 5.8,4 + 31.8a 



(12) 



Hron (1987) has used distances and radial velocities to 
young open clusters to investigate the rotation curve of the 
Galaxy. His results (valid for the range — 3 < i? — i? < 5 
kpc) are best fitted by A = 17.0 ± 1.5 km sec -1 kpc -1 
and a = —2.0 ± 0.6 km sec -1 kpc -2 . It has to be noted 
that Bochum 2 lies at <~14 kpc from the galactic center, 
therefore outside the limits of applicability for Hron's re- 
lations, but taking a = —1.4 km sec -1 kpc -2 (inside the 
quote errors) Eq.(12) gives a radial velocity of +68 km 
sec -1 , coincident with the observed value for Bochum 2 
(cf. Eq. 6). 
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